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ABSTRACT 


To study the nonlinear physics of inhomogeneous turbulent shear flow. 


the unaveraged Navier-Stokes equations are solved numerically. For initial 


conditions a three-dimensional cosine velocity fluctuation and a mean- 


velocity profile with a step are used. Although the initial conditions are 


nonrandom, the flow soon becomes turbulent. Concentrated turbulent energy 


develops near the plane where the mean velocity gradient is initially 


infinite. The terms in the one-point correlation equation for turbulent 


energy, including those for the diffusion and production of turbulence, are 


calculated; the diffusion terms tend to make the turbulence more homogeneous. 










The production, transfer between eddy sizes, and dissipation of turbulent 

1 o 

energy have been studied numerically. * In those references the turbu- 
lence is homogeneous, and no net spatial diffusion occurs. 

Here, the work is extended to an inhomogeneous developing shear layer so 
that diffusion, as well as the other turbulence processes, can be considered. 
The initial condition chosen for this purpose is, in dimensionless form 


Ui = ^ a" cos q"*x + ir6.j^V [sgn (x^ - n) + Vj 


where 

* 

^0 „ ^0 , n ^0 *n , ’^i 

u.= — u., V = — V, a.= — a., andx.=— . 

1 V 1 V 1 V 1 ^ ^0 

Note that the stars on dimensional quantities are omitted for corresponding 

dimensionless quantities. The subscripts can assume the values 1, 2, and 3. 

The quantity u^. is an instantaneous velocity component, V is a 

constant with the dimensions of a velocity, a^. is an initial velocity 

amplitude or Fourier component of the disturbance, q*"^ is an initial 
★ 

wavenumber vector, x. is a space coordinate, x^ is an initial 
characteristic length, v is the kinematic viscosity, «-jj is the Kronecker delta 
(equals 1 for i = j and 0 for i ^ j), and sgn x designates the sign 
of X. The quantities q^*x are dot products. The first term on the 


right side of Eq. (1) is the fluctuating part of the initial u-j, and the second 
term is the initial mean velocity. The latter is plotted against X 2 in the 
curve for t = 0 in Fig. 1. In order to satisfy the continuity condition 
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a?,? - 0 


For the present work 


= k(l,-2,l), a^ = k(l,-l,2) 


q} = q2= (1,1,1), qj = (1,-1,-1), 

where k is a quantity that fixes the Reynolds number. In addition to 
satisfying continuity, Eq. (4) gives 



ui = U2 = U3 = Uq 

at the initial time, where u. = u. - d. is a mean velocity component 

and the overbars indicate averaged values. Thus, Eqs. (1) and (4) give a 
particularly simple initial condition, as the components of the mean-square 
velocity fluctuation are initially equal. 
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The components in Eqs. (4) are the same as those in the earlier works ’ 
except that the signs of the a^ and q^ have been changed to make the 

initial u^^u^ negative. Thus, Uj^u^ does not have to change sign as a result 
of the dynamics of the flow, as it did in the earlier work, and the initial 
adjustment period is eliminated. If the adjustment period remained, much of 
the development of the shear layer would be distorted. 

To calculate the evolution of u^, the Navier-Stokes equations can be 
written in dimensionless form as 
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3U^ 3(U.U|^) 


~ 3^U. 

_ ^ + L_ 

3Xi 3X^3X^ 


where p is given by the Poisson equation 


j2p 3 




3X^3Xk 


t = -^ t*, and p = p* 


. * * 

The quantity t is the time, p is the pressure, and p is the 
density; the rest of the symbols are defined after Eq. (1). A repeated 
subscript in a term indicates a summation of terms. Equation (7) is 
obtained by taking the divergence of Eq. (6) and using the continuity 
equation, Eq. (2). 

In order to carry out the numerical solution of Eqs. (6) and (7) subject 
to the initial conditions of Eqs. (1) and (4) , a cubical grid with 32^ 
points and with faces at x^. = 0 and 2t is used. Modified periodicity 
is assumed for boundary conditions;^ that is. 


J 0 J 


^'^^x.=2w+b. ° ^'^^x.=b. 

W J V V 

for any bj. Equation (8) (not a tensor equation) is consistent with the 
initial condition given by Eq. (1) and is used to numerically calculate 
derivatives at the boundaries. 
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The spatial- and time-differencing schemes are essentially those used by 
Clark, et al.^ and in Refs. 1 and 2; that is, for the spatial derivatives, 
centered fourth-order difference expressions are used. For time- 
differencing a predictor-corrector method with a second-order predictor and 
a third-order corrector is used. The Poisson equation for the pressure is 
solved directly by a fast Fourier-transform method, which preserves 
continuity quite well. 

The calculated evolution of the mean velocity (U^ and U3 are 

zero) is plotted against in Fig. 1. The mean velocity is obtained by 

averaging over x^ and X3 at fixed values of X2. The shear 
layer grows (from essentially zero initial thickness) because of the presence 
of the turbulent and viscous shear stresses. The ratio of turbulent to vis- 
cous shear stress (averaged over x^ and X3 at the central plane 
X,, = tt) is plotted against dimensionless time in Fig. 2. Except at very 
early times the growth of the shear layer is almost completely dominated by 
the turbulent shear stress. 

Figure 3 shows the evolution of the instantaneous velocity component 
U2 and of the root-mean-square value of u^ (averaged over the central 
plane X2 = t^)- Although the initial conditions are nonrandom, the^^ 

evolution of U2 has a random appearance. On the other hand ^ evolves 
smoothly. These characteristics are representative of a turbulent flow. The 

quantity ^ increases monotonical ly at small times in contrast to the corre. 

spending curve in Ref. 1, where an initial adjustment period was present. As 
mentioned earlier (after Eqs. (4)), the initial adjustment period l^been 
eliminated here by changing the sign of the initial 03, so that u^u^ does 

not have to change sign as a result of the dynamics of the turbulence. The 
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decrease in 7 , near the end of the curve is caused by a decrease in mean 
velocity gradient, and thus of turbulence production, at large times (Fig. 1). 

As in the homogeneous easel, 5„,„_3cale fluctuations are generated in 
the inhomogeneous turbulence in fig. 3 by the interaction of the mean velocity 
with the turbulence. This can be seen by comparison of Fig. 3 with Fig. 1(a) 
in Refs. 1 and 2, where mean velocity gredients are absent. One might expect 
this, since it has been shown'’ that even for a general inhomogeneous tur- 
bulence, a term in the two-point spectral eguation for the turbulence can 
transfer energy between scales of motion as a result of the presence of mean 

gradients. 

A dimensionless plot of turbulent kinetic energy as a function of x; 
and time is given in Fig. 4. As for all of the averaged values. u,u„/2 is 
averaged over Xj and Xj for fixed values of x^. As time in- 
creases, an intense concentration of turbulent energy develops near the plane 
X 2 . ., where the mean velocity gradient is initially infinite. The tur- 
bulence is highly inhomogeneous. Inhomogeneity, in fact, seems to be the 
dominant characteristic of the turbulence generated in the shear layer. The 
indicated increase of turbulence with time is similar to that obtained ex- 

per imentally.^ 

Terms in the one-point correlation equation for the rate of change of 
the turbulent kinetic energy 
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are plotted for t = 0.000293 in Fig. 5. Equation (10) is constructed from 
Eq. (6) after letting Ui = Ui + and p = P + P- The terms that 

contribute most to the rate of change of v;;/2 are the production term 
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-lyjj dUj^/dX 2 , the pressure diffusion term -3pu2/ax2 and the 
kinetic-energy diffusion term -(1/2)3 vj^/3X2. The viscous diffusion 

term -3^C0j^/2)/3x^ and the dissipation term -3u,^/ 3 Xj^ 3 u^/ 3 X^^ are small 
in Fig. 5 . At early times, however, when the mean velocity gradient is 
large, the dissipation term is appreciable. 

The production term, whose form shows that turbulent energy is produced 
by work done on the Reynolds shear stress by the mean velocity gradient, is 
largest near the plane X 2 = ir, where the velocity gradient is initially 
infinite. The plots of the pressure and kinetic-energy diffusion terms show 
that those terms are negative near X 2 = if and positive away from that 
plane. Thus, they remove turbulent energy from the maximum. energy region 
and deposit it where the energy is smaller. Both diffusion terms therefore 
tend to make the turbulence more homogeneous. 

A comparison of the turbulence diffusion processes with the spectral 
transfer processes and the directional-transfer processes arising from 
the pressure-velocity correlations is instructive.^ The spectral-transfer 
processes remove energy from wave number (or eddy-size) regions where the 
energy is large and deposit it in regions of smaller energy. The directional- 
transfer processes remove energy from large-energy directional components 
and deposit it in a directional component (or components) where the energy 
is smaller. The turbulence-diffusion processes, as shown here, remove energy 
from regions of space where the energy is large and deposit it in regions 
smaller energy. The spectral-transfer, directional-transfer, and turbulence- 
diffusion processes tend, respectively, to make the turbulence more uniform 
in wave number space and more isotropic and homogeneous in physical space. 
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Although one might suppose that turbulence diffusion terms would always 
tend to make the turbulence more homogeneous, that supposition is not sup- 
ported by all experimental data. For instance, measurements of wall-bounded 
turbulence^ indicate that the pressure diffusion and the kinetic-energy 
diffusion terms transfer energy in opposite directions, although the net 
diffusion is from regions of high to regions of lower energy. On the other 
hand, measurements of turbulence in a free jet^ (closer to the case con- 
sidered here) seem to support the present findings. 
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Figure 3. - Calculated evolution of turbulent velocity 
fluctuations (normalized by Initial condition). Un- 
averaged fluctuations are calculated at center of 
numerical grid fX| * n). Root-niean-sguare fluc- 
tuations are averaged over Xj and Xj at central 

plane X 2 * m Uq' XqIv • 554, V ■ 2216 in 
Eq. (1). 
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